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Abstract 

We prove global existence and uniqueness of strong solutions to the loga- 
rithmic porous medium type equation with fractional diffusion 

^ + (-A) 1 / 2 log(l + n) = 0, 

posed for with nonnegative initial data in some function space of L logL 

type. The solutions are shown to become bounded and C°° smooth in (x,t) 
for all positive times. We also reformulate this equation as a transport equa- 
tion with nonlocal velocity and critical viscosity, a topic of current relevance. 
Interesting functional inequalities are involved. 
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1 Introduction 

In this paper we develop the basic existence, uniqueness and regularity theory for the 
problem 

' d t u+ (-A) 1 / 2 log(l + u) = 0, xeR,t>0, 

u(x, 0) = f{x) > 0, xeR. 

The equation in ( 11. li can be viewed as the limit m — > in the so-called fractional 
porous medium equation, 

(1.2) d t u + (-A) a/2 u m = 0, m > 0, < a < 2, 
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after a shift in the w-variable and a change in the time scale. The latter equation was 
treated in our papers [19], [20], where it was proved that it generates a contraction 
semigroup in L 1 (M Ar ) for any dimension N > 1, and that solutions become instanta- 
neously bounded and C a in space and time for data in L 1 ^) n L P (R N ) with p > 1 
larger than a critical value = N(l — m)/a. 

The difficulty we face here is that, according to those papers, the logarithmic dif- 
fusion is borderline for regularity questions when a = 1 and N = 1 for data in the 
natural space L 1 (R). This entails a very delicate critical-case analysis and a new type 
of regularity results: besides the expected result for / £ L 1 (IR) fl L P (M.) with p > 1, 
we obtain that solutions become immediately bounded when / belongs to an L logL 
space, almost L 1 (R) but not quite. This also offers some novelty when compared to 
the existing results for the standard porous medium equation, given by ( 11. 2ft with 
(j = 2, which are gathered in [31] and |32| . Actually, we go on to prove that the 
solutions are C°° in space and time, hence classical. 

Let us remark that the method proposed to tackle regularity has a more general 
scope. Actually, it can be applied to positive solutions of equations of the form 
d t u + (—A) a / 2 Lp(u) = posed in M. N under quite unrestrictive assumptions on the 
nonlinearity. We will study this issue in a forthcoming work. 

A further motivation for our study comes from the following connection: equation 
(II. ip can be transformed through a special nonlocal change of variables of the Back- 
hand type into the transport equation 

(1.3) d T v - H(v)d y v + d y H(v) = 0, y eR, r > 0, 

where H stands for a nonlocal operator which is a modification of the Hilbert trans- 
form. If instead of H we had the standard Hilbert transform H, using the identity 
d y H = (—A) 1 / 2 , valid when these operators are applied to regular functions, we would 

get 

(1.4) d T u - H(v) d y v + (-A) 1/2 v = 0, 

which is the transport equation with fractional diffusivity proposed by Cordoba, 
Cordoba and Fontelos in [10]. This is one of the several one-dimensional models 
considered in the last years to recast the main properties of the three-dimensional 
incompressible Euler equation and the two-dimensional quasigeostrophic equation, 
beginning with the work by Constantin, Lax and Majda [9]. 

Conveniently reformulated, our results for problem (11.11) produce existence and 
uniqueness of a classical global in time solution for equation (jl.3p for all initial data 
in L X (]R). This is a remarkable variation with respect to the results available for 
problem (II Ah : a global in time solution is only known to exist if the initial value 
belongs to the Sobolev space of fractional order if 1//2 (R), in which case it is in if 1 (IR) 
for almost every t > 0; see Dong [13j . From this regularity one might try to use the 
techniques of Kiselev, Nazarov and Shterenberg [15] to obtain further smoothness in 
space. 
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The application of the present approach to the transport equation f ll.4p is not im- 
mediate and needs further study. We believe that the connection between fractional 
diffusion and nonlocal transport problems is worth pursuing, since it may lead to the 
fruitful combination of very different techniques. 

As said above, the case a = N = 1 is critical in various aspects, in particular with 
respect to Sobolev embeddings. Thus, in the course of the proof of the smoothing 
effect in the mentioned LlogL space we need to use a critical fractional Trudinger 
type embedding due to Strichartz; see [27]. We generalize this embedding to other 
values of the exponents, a result in pure functional analysis that we hope could be of 
further application. 



2 Preliminaries and main results 

We recall that the nonlocal operator (—A) "/ 2 , o G (0, 2), is defined for any function 
g : M. N — > M. N in the Schwartz class through the Fourier transform, 

T{{-W* g ) (o = ier^)(o, 

or via the (hypersingular) Riesz potential, 

(2-1) (-^ /2 9(x) = C N , P.V. / dy, 

jr f ~~ y\ 

where Cn,o- is a normalization constant; see for example [16]. In our case, N = a = 1, 
the constant is — 1/tt, and we also have (— A) 1//2 = Hd x , where H denotes the 
Hilbert transform operator, defined trough 

Hf(x) = - P.V. / dy. 

If we multiply the equation in f 1 1.11) by a test function tp and "integrate by parts", 
we obtain 

(2.2) / ud^dxdt- / (-A) 1/4 log(l+ M )(-A) 1/ V^t = 0. 
Jo Jr Jo Jr 

This identity will be the basis of our definition of a weak solution. The integrals 
in (12. 2p make sense if u and log(l + u) belong to suitable spaces. The correct space 
for log(l + u)(-,t) is the homogeneous fractional Sobolev space ij 1 / 2 (IR), defined as 
the completion of C^°(1R) with the norm 



if 1 / 2 



L/2 =ll(-A) 1/4 < 



2- 



The Sobolev space if 1 / 2 (M) is then defined through the norm 

UsHMk + IK-A) 1 / 4 ^. 
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Definition 2.1 A function u is a weak L 1 -energy solution to problem (II. ip if: 

• u G C([0, oo) : L\R)) and log(l + u) G L 2 ((0, T) : HV 2 (R)) for every T > 0; 

• identity (12. 2p Zio/ds /or even/ y> G Cq(IR x (0, oo)); 

• m(-,0) = / almost everywhere. 

For the sake of brevity, we will denote the solutions obtained below according to 
this definition merely as weak solutions. We remark that this is not the only way 
of defining a solution to problem (II. ip . There are other possibilities, for instance 
entropy solutions [8], useful when dealing with equations involving convection terms. 

As for the initial data, our concept of solution only requires in principle / G L^(R). 
However, in order to prove existence we will ask / to belong to the slightly smaller 
LlogL-type space 

X = j/ > measurable : J {I + f) log(l + /) < oo j . 

Notice that L\(R) n L P (R) C X C L\(R) for any p > 1. This is to be compared with 
the result for the fractional porous medium equation (11.21) . where in the critical case 
a = N(l — m) we have required / G L 1 (IR Ar ) D L P (R N ) for some p > 1. 

The space X is natural for problem (jl.ip . Indeed, let ^(s) — (1 + s) log(l + s) — s. 
A function / belongs to X if and only if / G L^(R N ) and j R \&(/) < oo. On the other 
hand, after an integration by parts we formally obtain 

/ *(«(.,*))= / *(/)- r / K-A) i /4(i og( i +M ))| 2 < f 

and we conclude that the space X is preserved by the evolution. 
Notation. We will denote L x (f) := f R *(/) = / R ((1 + /) log(l + /)-/). 

Though we will be able to prove existence of a weak solution for any / G X, in order 
to prove uniqueness we will restrict ourselves to the smaller class of strong solutions. 

Definition 2.2 We say that a weak solution u to problem ( II. ip zs a strong solution 

if ^ue^oojxl). 

If it is a strong solution, then (— A) 1//2 log(l + it) is also an Lj^-function and the 
equation in ( 11. ip is satisfied a.e. 

Our first result shows that problem ( 11. ip is well posed in the class of strong solutions 
for initial data in X. 

Theorem 2.1 For every f G X there exists a unique strong solution to prob- 
lem fO]) . 
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Existence and uniqueness use an alternative formulation of problem (II. ip based in 
the Dirichlet to Neumann operator. Given a smooth bounded function g : R i— > R, 
we define its harmonic extension v = E(g) to the upper half-plane R 2 ^ as the unique 
smooth bounded solution to 

A Xt yV = 0, x 6 R, y > 0, 

f (x, 0) = g(x), x G R. 

Then, it turns out that — d y v(x,0) = (—A x ) 1 ' 2 g(x), where A^ is the Laplacian in 
all (x, y)-variables and A x acts only on the x- variables (in the sequel we will drop the 
subscripts when no confusion arises). The extension operator E can be defined by 
density in the space H^ 2 (R), and it is an isometry between this space and the space 
TL defined as the completion of C^°(R+) with the norm 



H 



JR 



1/2 

2 



Therefore, 



(2.3) / (-A) 1/4 (-A) 1 /^ = r f v E (0) • V EW. 
Jr Jo Jr 

We also have 

(2.4) / / VE(0) ■ VEty>) = / / Vt 7 -VE(V')- 



JO JO JR 

for any 7] eH such that Tr(r/) = 0; see |20] . 

Using this approach, problem ( II. ip can be written in an equivalent local form. If u 
is a solution, then u> = E(log(l + u)) solves 

( Aw = 0, (x,y) eR 2 + ,t> 0, 

(2.5) < d y w - dtP(w) = 0, x e R, y = 0, t > 0, = e 10 - 1. 

[ u; = log(l + /), xGR, j/ = 0, t = 0, 

Conversely, if we obtain a solution w to (12. 5p . then n = f3(w)\ Q is a solution to ( ll.ip . 
We next state the main properties of the solution obtained in the paper. 

Theorem 2.2 Let f G X . The unique strong solution u to problem (II. ip satisfies: 

(i) d t u G L 2 (R x (r, oo)) for all r > 0; 

(ii) Af-L 00 smoothing effect: 

(2.6) |K-,t)|U < CmaxirWpCcr 1 / 2 ^^)) 172 ),*- 3 / 4 !!/!!;/ 2 ^^)) 174 }; 

(iii) Lx{u(-,t)) and \\u(-,t)\\ p , 1 < p < oo, are non-increasing functions oft in 
(0,oo); 
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(iv) / u(x, t) dx = f(x) dx for every t > (conservation of mass) 
Jr Jr 



(v) u G C°°(R x (0,oo)); 



(vi) u(x, t) > for every x G R, £ > 0. 



PLAN OF THE PAPER. We will cover the existence and uniqueness theory in sections 
OandHl we borrow results and ideas from |2U] . Section is devoted to obtain some 
basic properties of the solutions. 

We then proceed with the smoothing effect, Section [61 first from L 1 fl L p , p > 1, to 
L°° and, then from X to L? . The proof entails a number of new ideas, in particular 
the use of a Trudinger inequality for fractional exponents. 

In Section [7] we perform a delicate regularity analysis to show that solutions are C°° 
smooth in space and time, and hence classical solutions of the equation. 

We next describe in Section [S] the transformation that passes from the equation 
in (II. ip to the nonlocal diffusion-transport model (II. 3p . and the results obtained for 
the latter. 

We finally include two appendixes. The first one is devoted to a generalization of 
the Nash- Trudinger type inequality used in the proof of the smoothing effect. In 
the second one we consider another tool used in that proof, an interesting calculus 
inequality. 

3 Uniqueness 

As mentioned in the introduction, in order to prove uniqueness we have to restrict 
the class of solutions under consideration. We will give two results in this direction: 
in the first one we restrict ourselves to weak solutions that satisfy u G L 2 (R x (0, T)) 
for all T > 0, and in the second to the class of strong solutions. 

Theorem 3.1 Problem (II. ip has at most one weak solution satisfying u G L 2 (R x 
(0,T)) for allT> 0. 

Proof. We adapt the classical uniqueness proof for porous medium equations due to 
Oleinik, Kalashnikov and Czou |18) . 

Let u and u be two weak solutions to problem ( II. ip . We subtract the weak formu- 
lations for u and u and take 




.og(l +u) — log(l + u))(x, s) ds, 



t > T, 



< t < T, 
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as a test function. Notice that, since the initial data of both solutions coincide, we 
do not need <£> to vanish at t — 0. After an integration in time we get 

T 

u — u)(x, t)(log(l + u) — log(l + u))(x, t) dxdt 







T \ 2 

(-A) 1/4 (log(l +u)- log(l + u))(x, s) ds) dx = 0. 







The condition m,!i£ L 2 (IRx (0, T)) ensures that the first integral is well defined. Since 
both integrands are nonnegative, they must be identically zero. Therefore, u — u. □ 

Remark. In particular, for / G L 1 (R) D L°°(R) there is at most a bounded weak 
solution. 

To prove uniqueness in the class of strong solutions we use the extension technique. 

Theorem 3.2 If u andu are strong solutions to problem (11. ip . for every < ti < t 2 
we have 

(3.1) / (u — u) + (x, t 2 ) dx < / (u — u) + (x, ti) dx. 



Proof. Let p be a smooth monotone approximation to the sign function such that 
< p < 1, and let j be such that f = p, j(0) = 0. Let ( e C^(R) be a cut-off 
function, < ( < 1, ((x) = 1 for \x\ < 1, ((x) = for \x\ > 2, and = ((x/R). 

Let z = log(l + u) - log(l + u). Using f l2T3|) and f )2^|) we get, for any < t x < t 2 , 
r diu-K± p{z)CR = _ f 2 f ( _ A)1 /4 Z ( _ A) i/ 4 (p(z)CR) 



dt 



t'2 POO p 

/ / VE(z)-V(p(E(z))E(( R )) 
ti Jo Jr 

t2 poo p 

/ / (p'(E(z))\VE(z)\ 2 E(( R ) + Vj(E(*)) -VE(Cfl)) 

< - r / (-a) i / 4 J (^)(-a)^ Cr < r / j^k-a) 1 / 2 ^! 

c f t2 f c(t 2 —ti) 

<- / \z\< max max{ u(-,t) i, u(-,t) i}. 

RJt! Jr r *6[ti,*a] 

where we have used that ((-A) 1 / 2 ^^)] = |(-A) 1/2 C(x/i?)|/ J R < c/i2, < j(z) < \z\, 
and the fact that log(l + u) < u for all u > 0. We end by letting R — > oo and p 
tend to the sign function. The case t\ — is obtained passing to the limit, using the 
^-continuity of u(-,t) at t — 0. □ 
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4 Existence of weak solutions 



The aim of this section is to construct a weak solution for any initial data in X. 
We will prove later, in section [5J that this solution, being strong, falls within the 
uniqueness class. 

Theorem 4.1 For every f G X there exists a weak solution u to problem (II. ip . 
This solution satisfies u > 0, 

(4.1) / / \(-A) 1 /Hog(l+ u)\ 2 dxdt<L x (f), 

Jo Jr 

and, iffeL°°(R), ||<,t)||oc < 



Proof. The construction of the solution uses several approximations. We refer to [20] 
for the details, where a similar calculation is made for the fractional porous medium 
equation (jl.2p . 

STEP 1. We first consider initial functions / G H L°°(R). We use the for- 

mulation of the problem in the extension to IR^_ version (12. 5p . By means of the 
Crandall-Liggett Theorem [12] we are reduced to deal with the elliptic related prob- 
lem 



(4.2) 



Aw = 0, x G M, y > 0, 

— d y w + (3{w) = g, x G M, y = 0, 



with 5- G L\{R) n L°°(R). Finally we substitute the half space ^ a ^alf ball 
5^ = {{x,y) : |x| 2 + y 2 < R 2 , y > 0}. We impose zero Dirichlet data on the "new 
part" of the boundary. Therefore we are led to study the problem 



Aw = in B%, 



(4.3) { w = on dB+ n {y > 0}, 

— d y w + /3(w) = g on := {|x| < R, y = 0}, 

with (7 G L°°{Dn) given. Minimizing the functional 



J{w) = \! \Ww\ 2 + [ (e w -(l+g)w) 
'b+ Jd r 
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in the admissible set A = {w G H l (B~l i ) : < (3(w) < \\g\\oo}, we obtain a unique 
solution it; = wr to problem (14. 3p . Moreover, if g\ and ^2 are two admissible data, 
then the corresponding weak solutions satisfy the L 1 -contraction property 

(/3(ioi(x,0)) - p(w 2 (x,0))) + dx < / (gi(x) - g 2 (x)) dx. 



D f 



STEP 2. The passage to the limit R — > 00 uses the monotonicity in R of the 
approximate solutions Wr. We obtain a function w^, = lim^oo Wr which is a weak 
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solution to problem (14. 2p . The above contractivity property also holds in the limit. 

Moreover, ||/3(io 00 (-, 0)) \\l°°{r) < \\g\\L°°(R), and > 0, since g > 0. 

Step 3. By the previous step, and using the Crandall-Liggett Theorem, we obtain 
the existence of a unique mild solution w to the evolution problem (12.51) . To prove 
that w is moreover a weak solution to problem (12. 5p . one needs to show that it lies in 
the right energy space. This is done using the same technique as in |19| . which yields 
the energy estimate 

III | VW(x, y, t)\ 2 dxdydt < L x (f) for every T > 0. 
Jo Jo Jr 

Hence the function u = /3(mJ(-,0)) is a weak solution to problem (II. ip . In addition, 
0))||i<x>( Rx ( 0;OO )) < ||/||l°°(r), and w > 0. In order to obtain estimate (14. ip we 
recall the isometry between i/ 1,/2 (IR) and "H. The Semigroup Theory also guarantees 
that the constructed solutions satisfy the /^-contraction property \\u(-, t)—u(-, t)\\i < 

11/ -/Hi- 

STEP 4. In this last step we consider general data / G X. Let {f k } C L^(lR)nL 00 (lR) 
be a sequence of functions converging to / in L 1 (M), and let {uk} be the sequence 
of the corresponding solutions. Thanks to the /^-contraction property we know that 
Uk{-,t) —> u{-,t) in L 1 (M) for all t > for some function u. Moreover, nonlinear 
Semigroup Theory guarantees that Uk — > u in C([0,oo) : [H]. On the other 

hand, using estimate ( I4.ip . we have log(l + Uk) G L 2 ((r, oo) : Z/ 1//2 (IR)) uniformly 
in k. Thus the limit u is a weak solution to problem ( II. ip for every t > r. The 
/^-contraction together with the /^-continuity allow to go down to r = 0. □ 



5 Strong solutions and energy estimates 

We still have to prove that the weak solutions that we have constructed are in fact 
strong. As a first step we consider the case of bounded weak solutions. The gen- 
eral case will follow by approximation as a consequence of the smoothing effect; see 
Section [6j 

Proposition 5.1 Let u be a bounded weak solution to problem ( II. ip . Then u is a 
strong solution and 

poo 1* 

(5.1) / / |^n| 2 rfx^<ct- 1 (l + ||n(-,t)|| 00 )L^(/), t > 0. 

Jt Jr 

Proof. In order to overcome the possible lack of regularity in time, we will work with 
the Steklov averages of functions g G L 1 1 oc (M x (0, oo)), defined as 

^ rt+h 

9 h (x,t) = - / g(x,s)ds. 
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A similar approach is used for instance by Benilan and Gariepy in [1] when dealing 
with evolution problems with standard Laplacians. The use of Steklov averages makes 
the process rather technical. The estimates are simpler to obtain when we assume 
regularity and work formally, and we invite the reader to do so. However, such 
regularity cannot be assumed at this stage of the theory. 

Almost everywhere we have 

M„ *\ - * ._ 9(x, t + h)- g(x, t) 



d t g h (x,t) = 6 h g(x J t) :~- 



h 



Let h > 0. Given any ip £ Cq°(M. x (0, oo)), we may take —S (p as a test function in 
the weak formulation. Then, using the "integration by parts" formula J °° J R <p 5 h u = 
— J °° J R u5~ h (p, we get that 

poo p poo p 

/ ip§ h udxdt = - (- A) 1 ' \\og(l + u)) h {- A) 1 dxdt. 

Jo Jr Jo Jr 

Taking ip = £d f (log(l + u)) h , where ( = ((t) £ C£°((0, oo)), this identity becomes 
(5.2) 

p oo p 1 poo p 

/ (d t u h d t (\og(l + u)) h dxdt =—z / C^|(-A) 1/4 (log(l + u))' l | dxdt 
Jo Jr * Jo Jr 

i poo p 

= - C |(-A) 1 /4(i g(i + u )ff dxdt. 

1 Jo Jr 

We now restrict ourselves to functions ( which are cut-off functions for the set [ti, t 2 }. 
To be more precise, we consider ip £ C°°(R) such that ifj' > 0, i/)(t) = for t < 1/2, 
ijj(t) = 1 for t > 1, and then define ((t) = ip(t/t{) — i])(t / (2t2)) . Then, using that 
C'(t) < t^ 1 maxip' , together with the inequality 5 h u 5 h log(l + u) > c(S h u) 2 , (with 
c = (1 + Halloo)" 1 ), we get 

c f 2 [ {5 h ufdxdt< ^~ I" I |(-A) 1/4 (log(l + M )) h | 2 dxdt. 
Jti Jr 2ti J J r 

The energy estimate (14. ip implies that the right-hand side is bounded for h small. 
Therefore there is a sequence h n — > + and a function g £ L 2 (M x (t, oo)) for all t > 
such that S hn u — > g weakly in L 2 (R x (t, oo)). It satisfies 

f [g 2 dxdt<ct^ 1 (l + \\u(-,t 1 )\\ 00 )[ [ |(-A) 1/4 log(l + u)\ 2 dxdt. 
Jti Jr Jo Jr 

On the other hand, 



ud t ip dxdt = — lim / I uS " ip dxdt 



oo 







lim / / 5 hn u<pdxdt 

h n ^0+ J Q 



gip dxdt, 
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which means that the distributional derivative d t u is in fact a function that coincides 
with g almost everywhere. □ 



We next prove that the ZAnorms do not increase with time. The main tool, used 
also later in the proof of the smoothing effect, Section 6, is the generalized Stroock- 
Varopoulos inequality |28| . [30j, 



(5.3) f A(z)(-A) 1/2 z> [ |(-A) 1/4 5( 

Jr Jr 



,2 



where A' = {B') 2 . An easy proof using the local realization of the half-Laplacian (in 
a more general setting) is given in [20l Lemma 5.2]. 

Proposition 5.2 Let u be a bounded weak solution to problem (11. ip . Then, for 
every < t\ < t2 we have 

L x (u(;t 2 ))<L x (u(;t 1 )), IK,* 2 )||p< \H-M)\\ P i l<P<oo. 

Proof. The first estimate is obtained directly multiplying the equation by log(l + u), 
as mentioned in Section [2j The cases p = 1 and p = oo in the second inequality 
follow from the elliptic estimates in Section HI For the rest of the cases, we put 
A(z) = u p ~\ z = log(l +u) in ([S3]). Since (-A) 1/2 z G L 2 (R) a.e. in t, if p > 3/2 we 
have A(z) G L 2 (IR). Assume this is the case. We then multiply the equation by A(z) 
and integrate in IR x (ti,t 2 ) to obtain 

- [ (u p (x,t 2 )-u p (x 1 t l ))dx<- f 2 [ \(-A) 1/4 G(u)(x,t)\ 2 dxdt < 0, 
V Jr v j J tl Jr 

where G{u) = B(z) = £ y/(p - 1) S p~ 2 /(1 + s) ds. 
For the case 1 < p < 3/2, we approximate the function A{z) by 



AJz) 



u p 1 for u > e, 
e p ~ 2 u for < u < e 



and then let e tend to zero. □ 
The L 1 -norm is not only non- increasing; it is conserved. 

Theorem 5.1 Let u be a strong solution to problem (jl.ip . For every t > we have 

u(x,t)dx= / f(x)dx. 



Proof. We take a nonnegative non-increasing cut-off function ijj(s) such that ip{ s ) = 1 
for < s < 1, ip(s) = for s > 2, and define 4>r{x) = ip{\x\/R). Observe that 
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\(-A) l / 2 (/) R (x)\ = J R- 1 |(-A) 1 / 2 ^(|x|/ J R)| < c/R. Multiplying the equation by <p R 
and integrating by parts, we obtain, for every t 2 > t\ > 0, 



u(x, t 2 ) — u(x, ti) ) 4>r(x) dx 



*2 

log(l + u)(x,t) (-A) 1/2 (f) R (x)dxdt 



ti 



< cR 1 max ||it(-,t)||i. 

te[h,t 2 ] 

In the last step we have used that log(l + u) < u. The result is then obtained just 
passing to the limit R — > oo. □ 

Weak bounded solutions turn out to have an energy which is well defined for all 
positive times. 

Proposition 5.3 Let u be a bounded weak solution to problem (II. ip . The energy 

E{t) = \ [ |(-A) 1/4 log(l + u)(x,t)\ 2 dx 

is a continuous function in (0, oo) which does not increase with time. Moreover, 
(5.4) E(t) < {2ty l L x (f) for every t > 0. 

Proof. Passing to the limit h — > in the identity ( 15. 2ft . we get 

l°° [ (\p^dxdt = - r [ C |(-A) 1/4 log(l + m)| 2 dxdt 
Jo Jr 1 + w 2 Jo J R 

for any test function ( £ C^°((0, oo)). This means that, as a distribution, E' coincides 
with the function — J R dx < 0. Since the latter belongs to L 1 (IR), we conclude 
that E £ VK 1,1 ((0, oo)), and therefore that it is a continuous function. Now we have 
L x {f) > 2 Jo E(s) ds > 2tE(t). □ 

In addition to the homogeneous Sobolev space HV 2 (R), the function log(l + u)(-, t) 
also belongs to the full space H 1 ' 2 { 



Proposition 5.4 Let u be a bounded weak solution to problem (II. ip . Then for every 
t > 

(5.5) || log(l + u)(;t)\\ H v> < t- 1/2 (LM)f 2 + ct-^\\f\\[ /2 (L x (f)) 1/A - 

Proof. Let w = log(l + u). We use interpolation and the Nash-Gagliardo-Nirenberg 
inequality (jA.6[) with N — 1, 7 = 1/2, q = 2, p = 1, to get 

\\<,t)h < M-M /4 M-M /4 < ciK-A) 1 /^^,*)!!^!!^.,*)!!;/ 2 . 

Next we use that log(l + u) < u, the energy estimate (15. 4p and the conservation of 
mass to conclude that 

(5.6) ||log(l + u)(-,t)|| 2 < C r^(L^ 



12 



□ 

To end this section, we improve the regularity of u and log(l + it), giving an L 2 - 
control of their gradients. 

Corollary 5.1 Let u be a bounded weak solution to problem (II. ip . Then u and 
log(l + u) belong to Lf oc ((0, oo) : H 1 ^)), and 

\d x \og(l + u)(x, S )\ 2 dxds < Ct"\l + \\u(;t)\\°o)L X (f), 

\d x u{x^)\ 2 dxds<cr l {l + \\u{-,t)\\oc?L x {f). 



oo 



Proof. It is clear from (JSSJ that log(l + u) E L 2 ((0,T) : L 2 (R)) for every T > 0. 
To estimate the gradient we just use (15. ip and the equation. Actually, 

/»oo /' /*00 /' 

/ / \d x log(l + u)(x,s)\ 2 dxds= / / \{-A) l / 2 log{l + u){x,s)\ 2 dxds 
Jt Jm. Jt Jr 

\d t u(x,s)\ 2 dxds < ct _1 (l + \\u(-,t)\\ OQ )L x (f). 
As to u, we just observe that d x u = (1 + u)d x log(l + u). □ 



6 Smoothing effect 

In Section H] we have constructed a weak solution of problem (jl.ip for general initial 
data / G A? by approximation with initial data in fl L°°(R). Our next aim is 

to prove that this solution becomes immediately bounded; in particular it is strong. 
Boundedness will follow from an estimate for bounded weak solutions, formula (]2.6p , 
which does not depend on the L°° norm of the datum, but only on Lx(f) (and time). 

The result will be obtained by combining L? — > L°° and X — > L 2 smoothing effects. 
The L? — \ L°° result is in fact a particular instance of a more general LP — > L°° result, 
valid for all p > 1. 

Theorem 6.1 Let u be a bounded weak solution, and let p > 1. There is a constant 
C > that depends only on p such that 

(6.1) !K,*)||oo < C max^/^H/llf p " 1) ,r 1/p ||/|| P }. 

We recall that the corresponding formula for the fractional PME with m > reads, 
in the case N = a = 1, 

(6.2) IK,<)lloo < cr l ^ m+ ^\\f\\f m+ ^ l \ 
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for every p > 1, cf. [2D]. Observe that when m = these exponents make sense for 
p > 1 but not for p — 1. It is also worth noticing that formula (16. ip can be obtained 
by formally putting in (16. 2p m = for u large and m = 1 for it small. 

Proof. The proof follows the same Moser iterative technique used in [20], but it 
is a little more involved. Let t > be fixed, and consider the sequence of times 
t k — (1 — 2~ k )t, = 2 fc p. We multiply the equation (recall that it is satisfied a.e. 
since u is a strong solution) by the test function 



Pk - 1 Pk 

and integrate in R x (t fc , tfc+i) (for p = P G (1> 3/2) we need an extra approximation 
argument, as in Proposition 15.21 to justify the computation). Using now the Stroock- 
Varopoulos inequality (15. 3p . we get 

Pfe+i 



Pfcbfe - 1) ' Pfc p/c(pfc 

>4 / ' / |(-A) 1/4 u Pfc/2 (x,r)| 2 rfxdr. 

Multiplying and dividing by ||ii(-, r)||£, for some r > 1, r > Pfc/2, using that the 
U norms do not increase in time, and applying the Nash-Gagliardo-Nirenberg type 
inequality flA.6[) with N — 1, 7 = 1/2, g = 2, we get 



(6.3) N-,fn-i)ii2£ < c2**" 1 ii«(-.**)ii;(n«(-,**)ii!;: + 



Let us denote [4 = max{||w(-, t k )\\ Ph , \\u(-,t k )\\ p Pk ^'' Pk }. Taking r = p k and r 



Pfc + 1 in (16. 3p we get that both t^+i) ||pfc+i an d IIM'; ^fc+i)llp^)+i are smaller than 
c2 k t~ 1 Uj^ k+1 , from where we obtain 

u k+1 < {c2 k r l fi^+^u k = (c2 k i p r l/p ) l/2k+1 u k . 

This recursive relation yields 

|K,t)||oo = lim U k < ct-V'Uo = ct-V'maxiWfW,, \\f\\^ )/p }. 

The final step is to get rid of the L p+1 -norm. Using Holder's inequality and the 
decay of the L p -norms we get 

IK,t)IL < c(t/2)^||/|| p max{l, \\u(; t/2) ||#}. 

If \\u(-,t/2)\\ 00 < 1, then ||w(-,t)|| 00 < c2 1 ' p t~ 1,p \\f\\ p and we are done. If, on the 
contrary, t/2)|| 00 > 1, we have 



<c(t/2)-^||/|L|| W (-,t/2)||^. 
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Since in this case, by the maximum principle, we have ||w( > ,r)|| 00 > 1 for every 
< r < t/2, we may iterate this estimate to get 

□ 

The above method does not allow to go down to p = 1. This drawback was already 
present in the PME case (both local and nonlocal, see [32] and [20]), where the limit 
exponent was p = max{l, (1 — m)N/a}. In the case N — a — 1, and putting m — 0, 
we get that the limit exponent should be p = 1, but it is not clear if solutions will 
become bounded when the initial datum only belongs to L 1 (M). Nevertheless, we 
may consider initial values in the slightly smaller space X. This is our next goal. 

Theorem 6.2 Let u be a bounded weak solution. There is a constant C > such 
that 

(6.4) f u 2 (x,t)dx < exp{c (r 1/2 (^(/)) V2 +r 1 / 4 ||/||; /2 (^(/)) 1/4 )} - 1. 

Proof. Fix any time t > and let w = log(i + w(-, t)). We know from Proposition 15.41 
that w G if 1/2 (M). Hence, using the Trudinger type inequality (1A.4I) . with N = 1 
and 7 = 1/2, we obtain 

I (> 2/c|Hl ^/ 2 _ ij < i. 

We now apply the calculus inequality (e w — l) 2 < (e fc — l)(e u,2//fc — 1) (see Lemma [B. II 
below for the proof), to get 

u 2 (x, t) dx < exp(c||w|| 2 ^i/ 2 ) — 1. 
We conclude using the energy estimate (15. 5p . □ 

To obtain the X — L°° smoothing effect we just have to combine Theorems 16.11 
andEja 



Proof of Theorem \2.2t (\\) . We first consider the case of initial data which are moreover 
bounded. The general case is dealt with by approximation. 

Using the LP — L°° estimate (16. ip with p = 2, and the X-L 2 estimate (16. 4p we get, 
first for t small, 

IK-,0IU < Cr l \\u{-,t /2)\\l < Cr'expiCt-^Lxif) 1 / 2 ), 

and then for t large 

IK-,i)lloc < cr^iH-^h < cr^\\f\\ l ' 2 L x {ff/\ 

Combining both estimates we obtain (12. 6p . □ 



15 



7 Regularity and positivity 



The solution that we have constructed in the previous sections is C°° for all positive 
times, and hence classical. This is the content of the present section. 



7.1 C ' a regularity 

The first and more difficult step is to prove that the solution u is C 1,Q for all a G (0, 1). 
Actually, given r > 0, u is uniformly C 1,a in Q T = R x (r, oo), denoted u G C^ ,a (Q T ) 
for short. 



Theorem 7.1 Let f & X. The strong solution to problem (11. ip satisfies u G 
Cu'^Qt) /or every < a < 1 and r > 0. 

Proof. STEP 1: u G (Qt) /or some a G (0, 1) and every r > 0. 

Once we know that u is bounded in the time interval t > r > 0, the result follows 
from the regularity results for problem ( 12.51) from Athanasopoulos-Caffarelli j3], since 
the nonlinearity (3(u) satisfies the non- degeneracy condition required in that paper. 

STEP 2: u G (Q T ) for every < a < 1 and every r > 0. 

To prove this, we will show that Holder regularity can be "doubled", following ideas 
from Caffarelli and Vasseur [6]; i.e., if u G C®(Q T ) for some a G (0,1/2), then 
u G C^ a (Q T ). The claimed regularity is then obtained repeating the argument a 
finite number of times. 

Let (xo,t ) G Q T be fixed, and denote u = u(xo,to). We write the equation in (11.11) 
as a fractional linear heat equation with a (nonlinear) source term, 

(7.1) d t u + M-A) 1/2 u = -(-A) 1/2 (log(l + u) - H- 

If we take fi = 1/(1 + Uq), the right-hand side of equation (17. ip can be written as 
-(- A) 1 / 2 F(n), where F(u) = log(//(l + u)) - fj,(u - u ) satisfies F(u ) = F'(u ) = 0. 
After a time shift, we may assume that u is uniformly C a and bounded down to 
t = 0. Recall now that the fundamental solution to the fractional heat equation 
d t u + (—A) 1 / 2 ?/ = is the Poisson kernel 

P(x,t) 



IT X 2 +t 2 



Taking a smooth approximation of P(x, fit) as a test function in the distributional 
version of (17. ip . and passing to the limit in the approximation we get that the solution 
u can be represented in the (mild solution) form 



u(x,t) = J P(x — x\, fj,t)f(xi) dxi 

(7.2) 



/ /(-A)^P 

JO JR 



(x — X\, fi(t — ti))F(u(xi, ti)) dxidtj 
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The first term in the right-hand side of (17.21) is regular, so we concentrate on the 
second one. 

We will use the notation y = (x, t) for the space-time variable, and also y = (x, /it) 
to accommodate the distortion in time created by the factor /x. We are thus led to 
study the regularity for the function 



(7.3) 



9(y) = A (y~ yx)X{t 1 <t}F{u(y 1 )) dy x 



where 



A(y) = A(x,t) = {-A)^P(x,t) 



1 x 2 -t 2 
n{x 2 + t 2 ) 2 ' 



Let us see first that the function g is well defined. To this aim we decompose M.+ as 
E p U Ep, where E p is the ellipse 

E = { Vl eR 2 + : - J7| < /£>}, 
with p small. Observing that 



1 



My-Vi)X{ti<t}dyi . 



x 2 - V 



dxdt = 0, 



we may write 



My - yi)x{ii<<} F ( M (?/i)) d yi 



< 



< c 



\A{y-y 1 )\\F{u{y 1 ))-F{u{y))\dy l 



dy x 



e \y - yi 



2-/3 



< C 



for some > 0, since both u and F are Holder continuous functions. On the other 
hand, using that F(u) is bounded we have 



Mv ~ yi)X{tx<t}F{u{yi)) dy A 



< c 



dx\dt\ 

V|a;i-a;|>l/(2/i) \ x ~ x l I 



< C. 



Now we will see that g(y) has the same regularity as F(u(y)). The key point is that 
if u is C a at yo, then F(u(y)) is C 2a at y$. Indeed, since F(uq) = F'(uq) = 0, and 
\F"(u)\ < c (recall that u > 0), we have 

\F(u{y))\<c\u(y)-u \ 2 <c\y-y \ 2a 

for every y G Ml. Moreover, the constants are independent of the point y . We ob- 
serve also that |^4(y)| < c|A(y)|, where c = c(/i). Since u is bounded and nonnegative, 
the constant c(/i) can be taken independent of /i. 
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Figure 1: Integration regions. 

Let y G be any point with \y — y \ — h. We have to prove that the difference 

(7-4) g(yo)-g(y)= [A(y -y 1 )x{t 1 <t }-A(y-y 1 )x{t 1 <t})F(u(y 1 ))dy 1 

Jr% v j 

is 0(h 2a ) for h small. In order to estimate the integral in ( 17.41) we decompose IR+ into 
four regions, depending on the sizes of | | and t\ — to, see Figure [TJ 

(i) The small ball = {\yi — yo\ < 4:h} C The difficulty in this region is the 

non-integrable singularity of A(y) &ty = 0. Integrability will be gained thanks to the 
regularity of F{u). We have, 

/ \A(y -y 1 )\\F(u(y 1 ))\dy 1 <c f < ch**. 

JB 4h J B 4h \yi ~~ Vol 

In order to estimate J B A(y — ^ 1 )F(m(|/i)) dyi, we consider as before the ellipse E c h, 
where c = c(/i) is chosen to have E c h C B 2 h, see Figure |2j 




Figure 2: Integration subregions in B^. 
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We get 

A{y ~ !7i)X{ti<i}-F(w(yi)) dyi 



A(V - yi)X{t 1<t} F(u( yi )) - F{u{y)) dy 



-a. 



,2a 



h 

+ F{u(y)) / A(y - yJxftKt} dyi, 

J B Ah -E ch 

V v ' 

h 

since the integral over the ellipse is again zero by symmetry. To estimate Ii we use 
the Mean Value Theorem applied to the function F to see that 

\F(u(y))-F(u( yi ))\ = \F'{6)\ \u{y)-u{ yi )\ < cm a x{\u(y) -u \, |u(j/i)-u |}|2/-2/ir, 
where 9 is some value between u(y) and u(y\). Therefore, 

\h\ < c [ , 1 J y - yx\ a (\yi - y \ a + \y- y \ a ) dyi < ctf 
JB ih \y - yi\ 

As to I 2 , since we are far from the singularity of A, 

|/ 2 | < c tf<* f ch ^. 

(ii) The narrow strip Sh = {\yi — 2/0 1 > 4/i, \t\ — t$\ < h}. In this region we have 
1 2/0 — 2/i I < § \y — yi I an d \xi — Xq\ > 3h. Therefore, 

\A(Vo -yi)X{ti<t } - A(y -yJxfaKtyWHuiViWdyi 

<f (l^(Fo - yi)\ + \A(y - \F(u( yi ))\ d Vl < r dyi 



s h v J Jsh Wo - Vi 

t0+h r dxxdt 



2-2a 



to-h J\x!-x \>3h \ X ~ %l\ 



(hi) The complement of the ball for large times, T^ = {\yi— yo\ > 4/i, t\ > t + h}. 
The integral in this region is 0, since here we have 

A (y - yik{ti<M = Mv - yi)x{t 1<t } = o. 



(iv) The complement of the ball B^ for small times, = {\y± — yo\ > Ah, t\ < 
to~h}. The required estimate is obtained here using the fact that we are integrating 
a difference of A's, so there will be some cancelation. Indeed, by the Mean Value 
Theorem, 

\A(y -y x ) - A(y-y x )\ < \y -y\ max{\d x A(Z)\, \d t A(C)\} < ch/\£\ 3 , 
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where £ = s(y — y ± ) + (1 — s)(y — y x ) for some s G (0, 1). On the other hand, since 
we are in D h , \y — y x \ < /i^ 2 |yo — < [£[ < c|£|, and we conclude that 

\MVo - yi)X{t 1<t0 } - A(y - yJxitKtyl < i °_ \ 3 - 
Therefore, assuming that a < 1/2, 

IMVo -y^X^Kto} ~ A(y -yjxit^tjl \F(u(y 1 ))\dy 1 

Dh r a 

<ch 

~ J Dh \yo-yi\ 3 ~ 2a ~ 

STEP 3: u G C^ a (Q T ) for every < a < 1 and every r > 0. 

We may assume, after a time shift, that r = 0. Let z = y — y$. The result will follow 
from an estimate of the quantity 

g(yo + z) - 2g(y Q ) + g(y - z) = / A{yQ,y,y 1 )F{u{y 1 )) dy x , 

Jm.1 

where 

A(y , y, yi) = A(y - !7i)x{ti<t} - 2A (Vo ~ !7i)x{ti<M + A(2y Q -y- y 1 )x{t 1 <2t -t}- 

As in the previous step, we consider separately the contributions to the integral of 
the four regions shown in Figure [TJ The contribution of the ball B^h is decomposed 
as the sum J x — 2J 2 + J 3 , where 



Ji= My - Vi)X{ti<t}F(u(yi)) dyx, 

■h= A(y - y 1 )x{t 1 <t }F{u(yi)) dy u 
Jb ah 

Ji= I A (rj - y 1 )x{t 1 <2t Q ~t}F(u(yi)) dy u rj = 2y - y. 

J Bah 

The integrals J\ and J2 were already estimated in Step 2. Since \fj — y \ = \y — y \, 
the integral J3 is estimated just in the same way as J\. 

The contribution of Sh is estimated in the same way as in Step 2, just using a rough 
estimate of the A's. The contribution of T/j is obviously 0. 

As for Dh, in this region we have, using Taylor's formula, 

ch 2 

\A(y ,y,yi)\ = \A(y - y ± ) - 2A(y - y x ) + A{2y - y - y x ) < — . 



12/0 — 2/i I 



Since 11 e C U °(R x (0, 00)) for all a G (0, 1), we obtain 



\A(y ,y, yi )\ \F{u{ yi ))\d yi < ch 2 I < ch 2a . 

D h JD h I2/0 - yi I 
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In summary we get 



g(y + z)- 2g(y ) + g(y - z) = 0(\z\ 2a ) 

for every a G (0, 1), \z\ < \yo\ (uniformly in y G K x (0, oo)). This estimate, together 
with the fact that g is bounded, allows to prove that (— A) a ^ 2 g(y) is bounded in Q T i 
for every a G (0, 2) and r' > 0. Indeed, if a G (0, 2a) and y G Qr', we have 

\(-&T /2 g(y)\ = 

f \z\ 2a , f dz 

^ C / U\2^ dZ + C / TT|2+^ r ^ C - 

>/{|*|<t'} H ■'{|*I>T / } l Z l 

Then, arguing in the same way as in |22l Proposition 2.9] (where the boundedness of 
the fractional Laplacian is assumed in the whole R 2 , not only in a half-plane), if we 
take a G (1/2, 1) and a G (1,2a), we obtain g G C 1,l3 (Q T i) for every G (0, a — 1), 
with uniform norm. We conclude that g G C\' a (Q T ) for every a G (0, 1), r > 0. □ 



7.2 C°° regularity 

Further regularity will now be a consequence of a result for linear equations with 
smooth coefficients which has independent interest. 

Theorem 7.2 Letv be a bounded weak solution to d%v + (— A) 1 ^ 2 (av + b) = 0, where 
the coefficients satisfy a, b G Cl> a (R x (0,oo)) n L°°(R x (0,oo)), a(x,t) > 0. // 
u G C"(R x (0, oo)) then v G C^ a (M x (r, oo)) /or ewen/ r > 0. 

Proof. Let (xo,to) G R+ be fixed and denote t> = v(xo,t Q ), a = a(xo,t Q ). Then v is 
a distributional solution to the inhomogeneous fractional heat equation 

d t v + a (-A) 1/2 v = (-A) 1/2 Fi + (-A) 1/2 F 2 , 

where 

Fx — -(a - a )(v - v ), F 2 = -b-v a. 
Reasoning like in the proof of Theorem 17. 1[ we are reduced to check that 

fi(x,t)= [ /(-A) 1/2 P(x-xi,a (t-ti))F 4 (xi,t 1 )rfx 1 c/ti, i = 1,2 
Jo Jr 

are C 1,a functions, with uniform norm. It is clear that f% inherits the regularity of 
F2; as to fx, we use the fact that the product (a — ao)(v — vq) is C 2a (or C x,2a ~ x if 
a > 1/2) whenever v is C a . □ 

Corollary 7.1 The strong solution to problem (II. ip belongs to C^(Q T ) for every 
t > 0. 



g(y + z) - 2g(y) + g(y - z) 



12+0- 



dz 
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Proof. The proof proceeds by induction. We know that u G C*> a (Q T ), a G (0, 1), 
t > 0. Assume that we have already shown that u G C^ ,a (Q T ) for some k > 1. Then, 
ffc = d^d2u, /3 + 7 = /c, satisfies an equation of the form <9 4 ?7fc + (— A) 1 / 2 (a/ c ffc + = 
0. Let us check that the coefficients satisfy the hypotheses of Theorem 17.21 On 
one hand, for all k > 1, df. — 1/(1 + u) is C h,a , hence C 1,a . It is also bounded, 
since u is nonnegative. On the other hand, as u G C^' a (Q T ) D L°°(Q T ), we obtain 
ffc G C°((5r) H L°°{Q T ). What is left is to verify that b k has the required regularity. 
In the case k = 1 we have bi = 0, and there is nothing to prove. When k = 2 we have 
three cases, 

(<%w) 2 (<9 x w) 2 , dtiid x u 



or 09 = -; — , or 



z (1+w) 2 ' z (1 + w) 2 ' (1 + w) 2 ' 

Since u G C 2 ' a (Q T ) n L°°(Q T ), we have clearly 6 2 G C^ a (Q r ) n L°°(Q r ). Applying 
Theorem Ol we obtain u 2 G C^ a {Q T >), t' > r. Hence u G C^*(Q T ,). 

The same reasoning works for every k G N. Indeed, the recursion formula for the 
coefficients b k has the form 

b k = dibk-x + v k -x dia , 

where i = x or i = t. We observe that b k is a polynomial in <9f <9/w, < f3' < /3, 
< 7' < 7, 1 < 0' + 7' < & — 1, with coefficients involving the powers (1 + u)~ m , 
< m < k. By the induction hypothesis, b k G C^ ,a (Q T ). As in the step k = 2 we 
conclude u G C^ +1 ' a (Q r /). □ 



7.3 Positivity 

Once the solution is regular and the equation is satisfied in the classical sense, we can 
use the Riesz representation (12. ip for the fractional Laplacian. Hence, at any point 
(xo,to) at which we have w(xo,to) = 0, we obtain 

Jr \Xq ~ S\ 2 

Since u is nonnegative, the right-hand side is nonnegative. Moreover, thanks to the 
conservation of mass, we know that the solution u is nontrivial if / ^ 0. Hence 
d t u(xo,to) is strictly positive. We have thus proved the following positivity result. 

Theorem 7.3 If f ^ 0, the solution to problem (II. ip is positive for all x G M 
and t > 0. 



8 A nonlocal transport equation 

We first recall that the half-Laplacian (—A) 1 / 2 can be written in terms of the Hilbert 
transform as (—A) 1 / 2 = Hd x = d x H. The latter equality holds provided that the 
operators are acting on a function belonging to some H^'^M) space, p > 1. 
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We now consider the change of variables (x, t, u) h- >■ (y, r, u) given by the Backlund 
type transform 



X 



y= (l + u(s,t))ds - c(t), r = t, v(y,r) = log(l + u(x,t)) 



o 



with c'(t) = ff (log(l + u))(0, t). We denote (y, r) = J(x,t). Notice that the Jacobian 

d(x,t) 



of the transformation J is d $ ,T } = 1 + ti ^ 0, since u > 0. Then we may write the 



inverse 

ry 



x = I e 

o 



-(^) da-c{r), 



with c'(r) = -#(log(l + u))(0,*)/(l + u(0,t)). 
We have 

d x y = l + u, d t y = -if (log(l + u)) = -if (v), 

where H(v) = H(v o J) o J -1 is the conjugate of the Hilbert transform if by the 
transformation J. Specifically, 

H(v(y, r)) = if (log(l + u(x, t))) = -P.V. / lo s(l + ^(<^)) d% , 

n Jr x — X 

= lp.v7 7F #i— rff 

With all this, equation (II. 1ft becomes 

(8.1) <9 r t> - if (u) <V' + ^ff (v) = 0, 

where i/ 6 1, r > 0. Since we assume u > we get u > 0. 
The L 1 norms of these two variables are related by 

u{x,t)dx= [ (1 - e- v{y > T) ) dy , 



and 

v(y,r) dy = / (1 + u(x, t)) log(l + u(x, t)) dx . 
Jr 

In particular, vq G L 1 (R) if and only if uq G X. On the other hand, 

d v v = d x u, d T v = — - — d t u + H (log(l + u))d x u. 

1 + u 

This allows to obtain regularity results for v from smoothness results for u. Finally, 
we have 



\(-A) 1/4 v(y,r)\ 2 dy = / l(-A) 1 / 4 log(l + u(x, t))\ 2 dx, 

Jr 

\d y v(y,r)\ 2 dy = [ \d x log(l + u(x,t))\ 2 (l + u(x,t))dx. 



Therefore, the results of the previous sections for (11.11) are translated to results 
for (18. ip as follows. 
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Theorem 8.1 Let v G L+(R). There exists a unique global in time classical solu- 
tion to equation (18. 1 P with initial value v . 

Theorem 8.2 Letvo G Li_(R). The classical solution v to equation ( 18. ip with initial 
value v satisfies: 

(i) L l -L°° smoothing effect: \\v{-, r) < C maxjr" 1 / 2 !!^!!^ 2 , t _3 / 4 ||?;o||i^ 4 for all 
t > 0; 

(ii) ||w(-,r)||i and t)||oo are non-increasing functions of r in (0, oo); 

(iii) (l — e~ v ( y ' T ^ dy = / (f — e~ v °^) dy for every r > (conservation law); 
Jr Jr 

(iv) v G C 1>a (]R x (0, oo)) for every < a < 1 ; 

(v) v(y, r) > for every y G M., r > 0; 

(vi) v G Lf oc ((0,oo) : if^R)). 

Appendix A: A Nash-Trudinger inequality 

In this appendix we contribute a new result that falls into the category of critical cases 
in embedding inequalities for spaces of functions with weak fractional derivatives. 

Theorem A.l Let <p G L P (R N ), 1 < p < oo, and assume that (-A) 7 / 2 G L q (R N ), 
< 7 < 1, q = N/^f. Put r = max{p, q], k = \p/r'~\, i.e., the least integer equal 
or larger than p/r', r' = rj(r — 1). There exists a constant a > such that if 
||0|| p +||(-Ar/ 2 0|| g <l then 

The particular case p = q < 2 was already proved by Strichartz in |27j . using 
estimates on Bessel potentials. Note that in this case k = 1, so that the integrand is 
just e a ^ P — 1. Our result covers all the possibilities for the parameters in the critical 
case. 

Before proceeding with the proof, we first review some related results and prelimi- 
naries for the reader's convenience. 

Sobolev SPACES OF integer ORDER. If 1 < p < N, Sobolev's embedding shows 
that W 1,P (R N ) is continuously embedded in L r (R N ) for all p < r < Np/(N - p). 
If p > N, then W l,p (R N ) is continuously embedded in L 00 ^^); even more, G 
C 0,1_iV ' p (R Ar ); the same happens for p = N = 1. The case p = N > 1 is critical and, 
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though W /1,7V (IR Ar ) L r (R N ) for every 1 < N < r < oo, it is easy to find examples 
of unbounded functions in H^ 1,Ar (M Ar ). However, if p — N > 1, then 



/ 



(A.1) 



N 



V 



N-2 

E 

i=o 



< 1 



/ 



for all in the unit ball of W 1,N (R N ), for some positive a independent of 0; see for 
example [I]. The proof of this result is based on the famous analogous estimate for 
the case of bounded domains due to Trudinger later improved by Moser, |17j . 

Fractional Sobolev spaces. Let l<g<oo,0<7<l. The homogeneous 
fractional Sobolev space W"?' q (R N ) is defined as the completion of C£°(M) with the 
norm 

||0lU 7 , q = ||(-A)^ 2 0|| g . 
The standard fractional Sobolev spaces are defined trough the complete norm 



■A) 



7/2, 



The well-known Hardy-Littlewood-Sobolev inequality 
Nq/(N- jq), then 

< C7||(-A)^/ 2 



[23] , states that, if g* 



q i 



for any 1 < q < iV/7 and < 7 < 1, and thus 



(A.2) 



W" r ' q (R N ) C 



Nq 



Bessel potential SPACES. The Bessel potential spaces are L^> q (R N ) = {/ = 
J 7 (0) : G L q }, 1 < q < 00, where J 7 (the Bessel potential of order 7 > 0) is defined 
in terms of its Fourier transform, 

XU)(0 = (i + l£lT 7/2 /(0- 

These spaces were introduced by Aronszajn and Smith [2] and Calderon [7], and have 
a natural norm, t \ \Li,q = || J_ 7 (/)||i«. The space L 1,q (R N ) is equivalent to the above 
defined fractional Sobolev space H /7 ' 9 (lR Ar ) for every < 7 < 1 and 1 < q < 00; see 
Stein (25. 

Critical Sobolev exponent. The inclusion (1A.2I) is not valid in the critical case 
q = N/j, as pointed out in [26]. However, if in addition we know that e L q (R ), 
then G L r (IR Ar ) for A^/7 < r < 00. That is, we have the inclusion 

(A.3) W J,N / y (R N ) C L r (R N ), for every A^/7 < r < 00. 

Indeed, in this situation belongs to the Bessel potential space L 1,N ^(R N ), and then 
the result follows from [27]. Notice that the case r = 00 is not included; see [26] . 
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To go beyond the //'-spaces, q < r < oo, in this critical case q = N/j, a careful 
estimate of the norms of the inclusion (IA.3I) using estimates of the Bessel potentials, 
allowed Strichartz [27] to prove the inequality 

(A.4) / f e °#l*^ -l)<l 

for some a > 0, valid for every such that 11011^7,^/7 < 1. That is, L 7 ' JV / 7 (R'' V ) is 
contained in the Orlicz space defined by the function in (IA.4I) . But this result is 
restricted to the range N/2 < 7 < N, unless the function has compact support. In 
our case of Sobolev spaces of fractional order < 7 < 1, this means that only N = 1 
can be considered, and then 1/2 < 7 < 1. 

BESOV SPACES. On the other hand, Peetre [2T] shows a restricted version of the 
previous inequality, in the spirit of (1A.1I) . valid for every < 7 < N, for functions in 
the Besov space A^ Nh (R N ); namely, there is a constant a > such that 

(A.5, / _£<!!!*«<!, 

k = \N/(N — 7)], for every such that ||0|| a jv/ 7 ,jv/ 7 < 1. 
The Besov spaces A^' 9 (IR Ar ) are defined through the norm 




A— = 

see [5]. It turns out that A^' 2 (M 7V ) = V4 /7,2 (lR Ar ). However, Besov spaces with p = q 
and Sobolev spaces are different whenever q ^ 2 [25 j. 

In the above-mentioned results, the control of both the function and its derivatives, 
or some quantity related to the derivatives, in the same LP space, yields a control 
in some Orlicz space. Our aim in Theorem IA.1I is to show an Orlicz-type estimate 
analogous to (|A.1|) and (1A.5|) starting from a control of the function and its derivatives 



in different LP spaces. 
We first obtain a generalization of the critical Sobolev-type embedding (1A.3I) . 

Proposition A.l Let p > 1 and < 7 < 1 we have 

L P (R N ) n W^' Nh (R N ) C L r (R N ) for every p < r < 00. 

Proof. It follows from the Nash-Gagliardo-Nirenberg type inequality 

\mi s <C(q n ,N)p\\(-Ay/ 2 ^M\ r P ~\ r = p + l-p/q, s = JV/(i\T - 7), 

valid for any function G L P (R N ) n W™(R N ), p > 1, q > 1, < 7 < 1, proved by 
the authors in [20]. Indeed, in the particular case q = A^/7 we have s = q', and thus 

(A.6) II0IIS:<C/||(-A) 7 / 2 0||«'||0||^. 
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□ 

Proof of Theorem \A.1\ As mentioned before, the particular case p = q < 2 was 
already proved in [27]. We will show how to treat the rest of the cases to get a 
complete analysis. 

Case p = q > 2. It was also covered in j27J under the additional restriction of 
asking to be compactly supported. For general functions some easy modification is 
needed. Indeed, for any function G L 7,q (M. N ), q = N/j, such that ||0||l7,9 < 1, the 
following estimate holds 

/ r \ l/H-W 

U\W<A[1 + - 



q' 

for every q < r < oo, where the constant A depends on N and q, but not on r, see 
|27| . We take then r = jq', j > q — 1 (which implies r > q), and obtain 

1^ —j^^ 7\ <oc 

j>q-i J ' i>-i J ' 

if we choose c > small enough. Finally, in order to have 1 in the right-hand side of 
( IA.1I) . we use that the function 



1 t jq ' 



m = * *-EV 



3=0 J 



satisfies F(Xt) < \ q 'F(t) for every t > 0, < A < I. 
CASE p < q. Using the Nash-Gagliardo-Nirenberg type inequality (1A.6I) we conclude 



that G L q (R N ), and thus G L 7,,? (]R Ar ). We apply then the previous case. 

CASE p > q. The key idea is that there is a value < fi < 7 such that G 
L At,7V//At (]R 7V ). Indeed, we can reach the exponent of integration p by lowering the order 
of differentiation. This follows from the Hardy-Littlewood-Sobolev inclusion (IA.2j) . 
which can be written as 

W^ N h?(R N ) c W^' Nhl (R N ) for every < ^ < 73. 

Hence, for the precise choice /i = N/p we obtain that (— A) M / 2 G ^(M^). We may 
now apply the case p = q with 7 replaced by fi = N/p. □ 

Remark. One is tempted to use the Nash-Gagliardo-Nirenberg inequality ( IA.6I) in 
order to estimate the sum in the development of the function in ( 1A.1I) . Unfortunately, 



the coefficient in flA.61) makes the sum divergent. 



Appendix B: A calculus inequality 

In the course of the proof of the smoothing effect we use a nice calculus inequality. 
Since it is not evident, we include a proof for the sake of completeness. 
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Lemma B.l For every x, a > we have 

(e ax - l) 2 < (e a - l){e ax2 - 1) . 

Proof. We develop the function f(x) = (e a — l)(e ax2 — 1) — (e ax — l) 2 in its Taylor 
series and rearrange the terms as follows: 

OO „ OO nl l, OO „ „ 

- a x a ( sr^ x a 



/w = E^E^- E 



n\ ^— ' fc! V / -— ' n! 

n=l fc=l n=l 

oo 

oo 

= V — - a n+k x n+k (x k ~ n - 1) . 
^ n\k\ 

By grouping the twin terms (n, fc) and (fc, n) we may restrict ourselves to the cases 
k > n and then 



oo oo 



f( x ) = J2 Yl ^ an+kxn+h ( xk ~ n + - 2 ) > °' 

n=l fc=n+l 

since the last factor is always positive for i^l and vanishes for x = 1. □ 
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